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I.  RESEARCH  WORK 

Research  in  the  report  period  continued  along  the  lines  indicated  in 
Progress  Report  I,  encompassing  the  following  areas: 

1.  The  analysis  of  the  kinetics  of  strongly  coupled  two-component  plasmas. 

2.  Study  of  the  ionization  equilibrium  in  a  strongly  coupled  plasma  medium. 

3.  Setting  up  a  scheme  appropriate  for  the  study  of  the  energy  levels  and  of 
the  ionization  potential  of  atoms  and  ions  in  a  strongly  coupled  plasma 
medium. 

1.  A  convenient  method  to  analyze  the  kinetics  of  strongly  coupled  plasmas 

is  by  way  of  partial  response  functions.  This  method  was  originally . suggested 

by  Singwi  and  Vashista1  and  then  further  analyzed  and  extended  by  Golden  and 

Kalman2.  However,  attempts  to  use  this  approach  for  concrete  calculations 

have  been  frustrated  by  the  apparent  complexity  and  unwieldiness  of  the 

formalism.  He  have  recently  succeeded  to  cast  the  formalism  in  a  much  simpler 

language  by  using  matrix  description  in  "species-space".  By  exploiting  this 

new  formalism  we  have  already  obtained  some  new  and  important  results: 

(a)  we  have  found  that  at  least  one  of  the  ion-ion  pair  correlation  functions 

in  a  multi-ion  system  exhibits  an  oscillatory  behavior  for  arbitrary  small 

coupling;  (b)  we  have  shown  that  if  quantum  and  proximity  corrections  are 

taken  into  account  in  the  interaction  potential  -  which  make  the  A-B 

interaction  different  from  the  A-A  interaction  (apart  from  scale),  -  the 

dielectric  function  exhibits  a  qualitatively  new  behavior;  (c)  by  examining 

the  structural  features  of  the  well-known  Totsuji-Ichlmaru  approximation  scheme 

for  multi-species  systems,  we  have  discovered  that  it  violates  an  important 

symmetry  requirement  and  therefore  its  applicability  to  multi-species  systems 

is  doubtful.  AIR  FORCEOfTICE  OF  SCIWTIF1C  RESEAT  MWfl 
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These  results  have  been  incorporated  in  a  paper  addressing  itself 
primarily  to  the  basic  formalism  (publication  (A));  further  papers  on  the 
results  listed  above  are  in  the  process  of  preparation. 

Another  aspect  of  our  research  is  based  on  an  approximation  scheme  (the 
Golden-Kalman  (GK)  approach) 3  which  combines  the  plasma  kinetic  equations  with 
nonlinear  fluctuation-dissipation  relations.  The  physical  transparency  and 
mathematical  tractability  of  the  resulting  chain  of  response  function 
equations,  which  has  already  provided  a  very  good  qualitative  description  of 
one-component  plasma  collective  modes  in  the  strong  coupling  regime^,  leads  us 
to  believe  that  the  GK  procedure  will  be  equally  effective  in  the  treatment  of 
the  corresponding  multi-component  plasma  problem.  Until  recently,  the  needed 
fluctuation-dissipation  relations  have  been  lacking. 

These  relations  have  now  been  established  by  us  for  both  binary  ionic 
mixtures  and  fully  ionized  ion-electron  plasmas.  Our  principal  frequency 
domain  result  links  a  single  three-point  dynamical  structure  function  to  a 
combination  of  nonlinear  partial  response  functions.  This  important  advance 
in  equilibrium  statistical  mechanics  takes  us  one  step  closer  to  acquiring 
information  about  transport  processes  in  the  more  realistic  bound  state 
plasmas,  for  it  will  now  be  possible  to  go  ahead  with  the  formulation  of  an 
approximation  scheme  leading  to  a  dielectric  response  function  description  of 
collective  modes  in  fully  ionised  multi-component  plasmas.  These  results  are 
given  in  publication  (B). 

Works  listed  under  (C)  and  (D)  deal  with  the  problem  of  plasma 
oscillations  in  strongly  coupled  systems.  Plasma  oscillations  are  expected  to 
reach  a  high  level  of  excitation  (~y)  for  strong  coupling.  The  ensuing 
dynamical  perturbation  on  the  bound  states  becomes  substantial  and  the 


knowledge  and  understanding  of  the  precise  plasmon  spectrum  under  these 
circumstances  is  essential. 

(C)  analyzes  the  relationship  between  perturbation  theoretic  and  the 
non-perturbative  GK  approach.  This  latter  aethod  was  applied  some  time  ago^ 
to  the  problem  of  plasmon  dispersion  and  gave  very  good  agreement  with 
molecular  dynamics  computer  simulations^  as  far  as  the  plasmon  dispersion  was 
concerned,  but  the  agreement  with  the  width  of  the  plasmon  peak  was  rather 
poor.  An  improvement  of  our  original  approach  described  in  (D)  now  has 
provided  a  remarkably  good  agreement  in  this  connection  as  well.  Including  the 
reproduction  of  the  non-monotonic  ^dependence  of  the  width.  To  the  best  of 
our  knowledge  this  is  the  first  successful  theoretic  description  of  this 
unexpected  behavior. 

2.  The  standard  method  for  determining  the  degree  of  ionization  in  a  plasma 
medium  is  through  the  use  of  the  Saha-equation.  However,  the  Saha-equation, 
based  on  a  model  of  a  gas  composed  of  independent  particles,  becomes  quite 
unreliable  at  high  densities.  The  main  effects  ignored  by  the  Saha  equation 
treatment  are  (i)  the  proximity  effect,  and  (il)  the  depression  of  the 
ionization  potential,  which  Itself  consists  of  the  depression  of  the  continuum 
and  of  the  lifting  of  the  bound  states  due  to  screening.  We  have  developed 
and  analyzed  a  model  which  shows  the  effect  of  the  latter.  The  influence  of 
the  medium  on  the  atom  is  described  through  the  replacement  of  the  intra- 
atomic  potential  by  a  screened  Debye  potential  where  the  screening  constant 
k  ■  (Aze^np)^2  depends  on  the  free  electron  density.  Increasing  values  of  k 
cause  the  disappearance  of  high-lying  bound  states,  leading  to  the  elimination 
of  the  last  bound  state  it  ica  ■  1.1  (i  ■  Bohr-radius)5.  Consequently, 
ionization  becomes  easier,  and  the  calculated  degree  of  ionization  is 
substantially  above  the  simple  Saha  result.  Moreover,  as  more  electrons 
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become  free,  the  ionization  potential  is  further  depressed  an  an  avalanche 
develops  generating  a  sudden  phase-transition-like  transition  from  the  weakly 
ionized  to  a  highly  ionized  state.  We  have  used  the  recent  analytic 
approximation  given  by  Green^  to  describe  the  bound  states  in  a  screened 
potential.  The  Saha  equation  was  modified  by  using  the  Larkin  correction7  for 
the  partition  function.  The  resulting  rather  complex  self-consistency 
condition  for  the  degre  of  ionization  was  solved  numerically.  Detailed 
couqmter  work  has  been  performed  recently  with  parameters  pertaining  to 
U-pla8mas  in  a  wide  range  of  densities,  up  to  the  degeneracy  limit.  We  have 
found  the  phase-diagram  in  the  density-temperature  plane,  with  a  critical 
point  in  the  vicinity  of  Tc  *  3  eV,  n,.  -  3  x  10^3  cm~3.  These  results  are 
described  in  publication  (E).  Further  works  on  ions  with  higher  Z  values,  and 
on  the  effect  of  the  electron  degeneracy  are  in  progress. 

3.  There  are  a  number  of  methods  available  for  the  study  of  the  dynamics  of 
an  atom  or  ion  immersed  in  an  ionized  medium7.  The  Thomas-Fermi  (TF)  method 
combines  relative  simplicity  with  excellent  qualitative  descriptive 
capability.  While  the  TF  scheme  and  its  refinement  have  been  employed  to  the 
problem  on  hand,  the  approaches  used  are  appropriate  for  weakly  coupled 
situations  and  not  for  the  case  of  strong  coupling.  The  scheme  we  have  worked 
out  is  intended  to  include  specific  strong  coupling  effects  and  is  summarized 
below. 

-  The  system  is  considered  to  consist,  for  the  purpose  of  the  TF  equation, 
of  three  species: 

(i)  Ni  classical  nuclei  of  charge  Z, 

(li)  Ne  -  N^X  hot,  classical,  "free”  electrons,  and 
(ill)  Njj  *  N^Y  (Y"Z-X)  degenerate  bound  electrons. 


-  The  potential  around  a  chosen  nucleus  is  determined  by  the  Poisson- 
equation  with  the  and  gn  providing  the  density  distribution  of  the 
surrounding  particles. 

-  For  the  purpose  of  the  calculation  of  the  correlation  functions  the 
system  consists  of 

(i)  ions  of  charge  X  and 

(ii)  Ne  hot  classical,  "free”  electrons.  Interacting  through  the 

effective  potentials  *ee>>  *ii  and  *ei»  such  that  z  >  zeff  >  x* 

-  Zef t ^  and  Zeff,  e  are  to  determined  by  the  integrated  charge 
densities  of  the  bound  electrons  only. 

-  The  correlation  functions  gee  and  ge^  are  now  to  be  determined  by  using 
one  of  the  strongly  coupled  static  plasma  schemes-elther  the  mean  field 
theories,  or  HNC. 

-  Finally  X  is  to  be  determined  by  calculating  the  free  energy  of  the 
system  and  minimising  it  with  respect  to  X. 

The  salient  feature  of  this  scheme  is  that  it  applies  the  philosophy 
of  the  TF-DH  (Debye-Huckel)  method,  but  replaces  the  OH  distribution  with 
the  correlation  functions  appropriate  for  strong  coupling.  In  turn,  the 
correlation  functions  depend  on  the  effective  potentials. 

We  have  prepared  a  survey  of  the  earlier  TF  and  TF-DH  works,  starting 
with  the  early  works  of  Feynman,  Metropolis  and  Teller**  and  Latter^  through 
the  more  recent  works  of  Rossnyai  and  Alder10  and  More  and  Skupsky11,  which  we 
plan  to  make  available  as  a  Technical  Report.  The  superiority  of  our  scheme 
as  outlined,  over  the  existing  approaches  is  quite  convincing. 

There  are  a  number  of  theoretical  problems  which  have  to  be  clarified 
before  detailed  numerical  work  can  be  performed: 


(i)  Che  proper  definition  of  the  chemical  potential  and  chemical 
equilibrium  for  the  two-electron  mixture; 

(ii)  the  behavior  of  the  mean  field  theories  for  a  multi-species 
system  (see  above); 

(iii)  the  correct  combination  of  the  mean  field  theory  approach  with  the 
TF  scheme. 

These  questions  are  presently  being  investigated. 
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II.  FUTURE  PLANS 


We  will  continue  to  examine  the  three  major  problem  areas  listed  above, 
with  increased  eaq>hasis,  however  on  topics  (2)  and  (3),  since  the  main 
objective  of  (1),  namely  the  extension  of  the  ocp  formalism  to  the  case  of  the 
multi-species  systems  has  been  accomplished.  As  to  (2),  we  will  include 
quantum  (degeneracy)  effects  in  the  deformation  of  the  phase  diagram  and  will 
generalize  for  arbitrary  Z  values.  Furthermore,  we  will  attempt  to  analyze 
the  depression  of  the  continuum  (as  contrasted  to  the  lifting  of  the  bound 
states)  -  which  is  a  rather  difficult  problem.  As  to  (3),  we  hope  to  be  able 
to  clarify  the  pending  theoretical  problems  and  to  start  actual  numerical 
calculations. 

A  further  topic  to  be  investigated  is  the  effect  of  plasma  oscillations 
on  the  bound  states.  We  will  analyze  the  level  of  oscillations  from  the 
available  computer  simulation  data  and  will  examine  the  possibility  of  the 
merging  of  the  continuum  with  the  high-lying  bound  states. 
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